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Abstract 

Relational particle models are of value in the absolute versus relative motion debate. They are also analogous to the 
dynamical formulation of general relativity, and as such are useful for investigating conceptual strategies proposed for 
resolving the problem of time in quantum general relativity. Moreover, to date there are few explicit examples of these 
at the quantum level. In this paper I exploit recent geometrical and classical dynamics work to provide such a study 
based on reduced quantization in the case of pure shape (no scale) in 2-d for 3 particles (triangleland) with multiple 
harmonic oscillator type potentials. I explore solutions for these making use of exact, asymptotic, perturbative and 
numerical methods. An analogy to the mathematics of the linear rigid rotor in a background electric field is useful 
throughout. I argue that further relational models are accessible by the methods used in this paper, and for specific 
uses of the models covered by this paper in the investigation of the problem of time (and other conceptual and technical 
issues) in quantum general relativity. 
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1 Introduction 



In Euclidean relational particle mechanics [Tl[2l[3l|4l[5l[6l[71[8l[9], only relative times, relative angles and relative separations 
are meaningful. While, in similarity relational particle mechanics [lOl lU [6l [7l [TTl [12] , only relative times, relative angles 
and ratios of relative separations are meaningful. These mechanics are valuable models as regards the absolute versus 
relative motion debate [13j . It is then of interest what structure one gets when one quantizes such theories. There has not 
been much work on this to date as regards nontrivial explicit examples [ |14| concerns formal quantum constraints, [15" is 
a toy of geometrical quantization, [161 [3] is a semiclassical toy (but with only explicit examples for Euclidean relational 
particle mechanics in 1-d) and some simple solutions of the Dirac quantization scheme [5]]. This paper goes further 
than these works in being the first quantum treatment of explicit nontrivial examples of similarity relational particle 
mechanics, while \1T goes further, rather, by considering more complicated explicit examples of Euclidean relational 
particle mechanics than those considered before. We get there using our recent understanding of reduced configuration 
spaces [11] as a path to quantization. 

In investigating conceptual strategies for the Problem of Time ITS', 'T9', "SOl [HI [52 [2l [25l [1 [26l [27] in Quantum Gravity, 
relational particle mechanics are useful analogues of GR [531 [2 [24, 3, 26j[S|; I view this as major long-term motivation 
for the current series of papers. This notorious problem occurs because 'time' takes a different meaning in each of GR 
and ordinary Quantum Theory. This incompatibility underscores a number of problems with trying to replace these two 
branches with a single framework in situations in which the premises of both apply, namely in black holes and in the very 
early universe. One facet of the Problem of Time that shows up in attempting canonical quantization is that the lack of 
linear momentum dependence of the GR Hamiltonian constraint leads to a frozen (i.e. timeless, or stationary) quantum 
equation for the universe arises therein: the quantum counterpart of (I-H) is the Wheeler-DeWitt equation 

EE |VA1AA^''''^^| - eRic(X)| 'f - Vh{Ric{h) - 2A}* + 7i,„,ttc,* = , (1) 

where 5* is the wavefunction of the universe. 

Before further consideration of the Problem of Time, a detour is first required about other features of this equation 
and of the toy models of it that are the subject of this paper. The inverted commas indicate that the Wheeler-DeWitt 
equation has, in addition to the Problem of Time, various technical problems, including 

A) regularization problems - not at all straightforward for an equation for a theory of an infinite number of degrees of 
freedom in the absense of background structure, while the mathematical meaningfulness of functional differential equations 
is open to question. I emphasize that this is not an issue in this paper's toy models as these are for a finite number of 
degrees of freedom, similarly to the situation in minisuperspace quantum cosmology. 

B) There are operator-ordering issues, which this paper's toy models do exhibit an analogue of. To view the analogy, let Qa 
be general coordinates (spanning spatial indices and either particle labels or field species along with spatial dependence) 
with a corresponding configuration space metric A^ab with inverse Af^^ and determinant A4 (perhaps merely at the 
formal level) . Then let be a partial derivative for finite systems or (perhaps merely formally) a functional derivative 
for infinite systems. Then the Laplacian ordering for the classical combination of configurations and their momenta 
J\f''^{Q^)VAVB is 

D^= (VtwAA^^-^I , (2) 



VA^VQa I VQe 

which has the desirable property of being independent of coordinate choice on the configuration space ^28J . This property 
is not, however, unique to this ordering: one can include a Ricci scalar curvature term so as to have — ^Ric(Al) 
[551 [551 [3D] . There is then a unique conformally- invariant choice [301 [3T| among these orderings. The conformal invariance 
here corresponds to retaining at the quantum level the banal conformal invariance that is obvious and natural in the 
relational action at the classical level [H [31] . This ordering is 

V'= f^j^'^B^X- J' -^ mciM) , (3) 

VMVQaX VQb J 4{/e-l} ^ ^' 

where k is the configuration space dimension. Furthermore, for this to be conformal, it is required that ^' itself transforms 
in general tensorially under the conformal transformation [32 j Fl 

<If — >^i^ = n^^ . (4) 

Moreover, this paper focusses on a model with a 2-d configuration space, for which the conformal £, — {k — 2}/4{fc — 1} 
collapses to zero, so that Laplacian ordering and conformally invariant wavefunctions suffice (but almost all other relational 
particle mechanics models, such as |17j . have configuration space dimension > 3 for which this subtlety is required). Finally, 
if one sends H^* — E^f to Hn^E'n = En\l/n = {E/il'^l^'n, one has now an eigenvalue problem with a weight function 



^This paper's models, like minisuperspace, only require the finite configuration space dimension k case of these equations. 
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which then appears in the inner product: 



{■^i)n*{-^2hn-^VM^d''x . (5) 

This inner product additionally succeeds in being banal conformally invariant, being equal to 

/ ^'l*^]^*2^^^^^~^V7W^^''d'^a;= / *i**2%/A^d'^a; (6) 

in the 'original or mechanically natural' representation in which E comes with the trivial weight function, 1. I should 
caution however that, while this conformal choice of ordering is a nice choice which I am suggesting lies on a more solid 
principle than how it is usually presented (fully explained in [31]), it might nevertheless eventually be found to clash with 
other requirements such as existence and suitable behaviour of crucial operators. 

C) Next, and talking only in the context of finite models such as in this paper rather than for full geometrodynam- 
ics for which the mathematical machinery is lacking^ in general Hq is not self-adjoint with respect to n{ \ )a, while 
the mechanically-natural H is, in a simple sense, with respect to ( | ). I.e. in the sense that / \/7Ud''x**L>2^ ^ 
/ VMd''x{D^'^*}^'^ + boundary terms, which amounts to self-adjointness if the boundary terms can be arranged to 
be zero (which is definitely not a problem in this paper as there are no boundaries) and in other cases involves such as 
suitable fall-off conditions on 5*. This is not shared by the fi-inner product as that has an extra factor of fi^^, which in 
general interferes with the corresponding move by the product rule (VAT does not interfere thus above, since the Lapla- 
cian is built out of derivatives that are covariant with respect to the metric A^ab-) However, solving Hn^o = Eo^'a is 
equivalent to solving H^E' = E^I*, so the banal conformal transformation might at this level be viewed as a sometimes- useful 
computational aid, with the answer then being placed in the mechanically-natural representation for further physical 
interpretation. This is not an issue in this paper as O is but 1/4 in my spherical calculations, thus not presenting any 
product-rule obstacles to self-adjointness. 

Now, returning to the Problem of Time, while many conceptual strategies have been put forward to resolve it, there is 
a long history of proposed resolutions not standing up to detailed examination 23 , so that this remains an open problem 
for GR. Some of these strategies are as follows. 

1) Perhaps within GR at the classical level there is a fundamental hidden time. E.g., one could seek for such a time 
by canonically transforming the geometrodynamical variables to new variables among which an explicit and genuinely 
time-like time variable is isolated out. One candidate time of this form is the York time, [35112011^. This is proportional 
to hp^iyTT'^'^ / ^/h so it is a 'dilational object'. 

2) Perhaps instead there is no fundamental time in Quantum GR but a notion of time emerges in the quantum regime e.g. 
in regions of the universe that behave semiclassically [26l[36]. In situations in which the Born-Oppenheimer approximation 
\1/ = V'(Ha')Ix(Ha', La")) for Ha' 'heavy, slow' and La" 'light and fast' degrees of freedom and the WKB approximation 
^/'(Ha') = e*'^'^*'-'/^ are applicable, an emergent time drops out of the WDE 19, 36, 26l[37]. For (e.g. using hf^i, as 'H' 
and the matter as 'L', and requiring that \x) depends nontrivially on H so that the QM is rendered nonseparable in H, L 
variables), 

^2^p.pa ^ ^ contains /i^AA^'-"'^-!^!!^ = (7) 

(by identifying F « \N /M where W = W(Ha') is Hamilton's principal function and M is a generic H-mass, and then using 
the Hamilton-Jacobi relation for the momentum and 
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Then by the momentum-velocity relation (1. 10), this expression contains 



(8) 



by the chain-rule in reverse so that one has a TDSE for the local L degrees of freedom with respect to a time standard 
that is (approximately) provided by the background H degrees of freedom. An issue here is that (semi)classical conditions 
need not always occur ~ guarantee of a classical 'large' as in the Copenhagen Interpretation of QM has been cast aside 
in considering the universe as a whole, and will then by no means be recovered in all possible situations. The above 
semiclassical approach is additionally a useful framework for discussing the origin [36' of galaxies and cosmic microwave 
background perturbations within the semiclassical scheme, for which one needs to study spatially-located fast light degrees 
of freedom that arc coupled to global slow heavy degees of freedom such as the size of the universe. 



■^The Loop Quantum Gravity approach 1331 is advantageous at this point in being equipped with a Hilbert space structure, as well as being 
better-suited as regards A) above. Its passage to Ashtekar variables renders operator-ordering issues there different to those of geometrodynamics 
and the toy model of this paper. Nevertheless, discussion of analogies with the geometrodynamical Wheeler-DeWitt equation as regards the 
Problem of Time remain reasonable and commonplace in the literature, e.g. 1231 1341 1261 make use of such. 
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3) There are also timeless records strategies [HI [Ml HI Ull ISl 1321110] • Here, the primary objects are records - information- 
containing subconfigurations of a single instant that are localized in both space and configuration space. One would then 
seek to construct a semblance of dynamics or history from the correlations between such records. QM probability density 
functions on shape space are here of interest as regards whether one can substantiate Barbour's conjectures [24l [3] about 
a present populated with time-capsules. 

4) There are also approaches in which it is the histories that are primary [351 HI] ■ 

The currently intended applications of this paper's toy models are to approaches 2), 3), 4) and combinations thereof. 
However, for completeness and due to current interest, I also outline the following additional family of strategies. 

5) Distinct timeless approach involve evolving constants of the motion (a Heisenberg rather than Schrodinger type ap- 
proach) or partial observables [34 , which is used in Loop Quantum Gravity's Master Constraint program 33J . 

Quantum GR being technically difficult, toy models such as that being developed in the current paper have been useful 
toward developing strategies such as the above. Relational particle mechanics are useful such, due to geometrodynamics 
being parallely formulable in relational terms (see [42l [43] and [42l |44] for its robustness to the inclusion of matter, 
or for a summary. Sec 1.2.3-4). Specific Problem of Time applications and extensions of this analogy then feature in 
[TBI l23l [21 [24l m m [SI [6l [37l [45l [46l [47l [40] ) . In particular, both relational particle mechanics and geometrodynamics have 
a constraint that depends quadratically but not linearly on the momenta, which feature underlies the frozen formalism 
aspect of the Problem of Time, and both have further nontrivial constraints that are linear in the momenta, which cause 
many of the complications with strategies proposed to resolve the Problem of Time. Relational particle mechanics have 
so far been useful toy models in possessing an analogue of the above-mentioned dilational York time [S] [371 llS]i for 
the semiclassical emergent time approach [371 ESI 113 ^^'^ for the timeless records theory approach [51 [211 [3] HO] . See the 
Conclusions of this paper and [17^ for brief discussion of useful applications of (extensions of) these papers' models to these 
issues. In relational particle mechanics models, by their constraints and the subsequent appearance of more complicated 
reduced configuration spaces on which these models' partial observables would live, the partial observables approach to 
these examples of nonrelativistic mechanics would be rather distinct from that for the mechanics models considered in 
|34| . and may lead to a broader view than therein on what features distinguish nonrelativistic mechanics theories from 
relativistic ones. 

Both the absolute versus relative motion debate and the study of conceptual strategies suggested toward resolving the 
Problem of Time in Quantum Gravity benefit from study of explicit examples of quantum relational particle mechanics. I 
provide such in this paper (similarity relational particle mechanics) and [17 ] (Euclidean relational particle mechanics) , by 
applying my recent understanding of reduced configuration spaces |llj to carry out quantization. [The above-mentioned 
QM investigations of relational particle mechanics [TH [15l [TBI [5] , and the recent preprint [48] (path integral method), 
carry out their quantization by other means.] In this paper I provide Schrodinger equations for scalefree N-stop metroland 
(relational mechanics of N particles in 1-d) and scalefree N-a-gonland (relational mechanics of N particles in 2-d), concen- 
trating on the latter 's triangleland case (3 particles in 2-d). For this I discuss how relative angle independent potentials 
are a substantial simplification (in close analogy with how central potentials are in ordinary QM) , and provide some simple 
solutions. I begin with relative angle independent multiple harmonic oscillator type potentials, for which I additionally use 
asymptotic, perturbative and numerical methods. It is important for this work that I identified my problem to have the 
same mathematics as the Stark effect for the linear rigid rotor. In Sec 3 I use rotated coordinates/normal modes/adapted 
bases to remove the relative angle independence restriction: any multiple harmonic oscillator for scalefree triangleland 
admits coordinates in which it takes the form of Sec 2's special case, which amounts to the 'electric field' being in an 
arbitrary direction rather than in the simplifying adapted basis in which it points along the Z-axis (of the embedding 
of the configuration space sphere) . This is of importance as regards setting up models of the semiclassical approach and 
records theory. My Conclusion (Sec 4) includes an outline of several further examples to which techniques of this paper 
can be applied and which are also useful as regards modelling the Problem of Time and Quantum Cosmology. 



2 Quantum Similarity Relational Particle Mechanics 

2.1 Time-independent Schrodinger equations for Scalefree N-stop metroland and N-a- 
gonland 

For scalefree N-stop metroland, the configuration space is S"~^ and the conformal-ordered time-independent Schrodinger 
equation is (via App A.l and the notation of Sec 1.3.2 and with E being redefined to absorb the conformal contribution, 
which is possible as spheres are spaces of constant curvature) : 

"e^.-^l + v^^E^' . (10) 



2 sin"-i-^"e^ n^;} sin^Gp 56, I 96, 

For scalefree N-a-gonland (N > 3), the configuration space is CP"^^ and the conformal-ordered time- independent Schrodinger 
equation is (via App A.l and likewise being able to absorb the conformal contribution as complex projective spaces are 
Einstein and hence of constant Ricci-scalar curvature |11| ) 
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- - 9* 



V* = E^- . (11) 



2.2 Time-independent Schrodinger equation for scalefree triangleland 

In this paper we mostly consider the triangleland case for which there is the additional 'accident' CP^ ~ which permits 
use of both spherical type and complex projective type variables and carries stronger guarantees of good mathematical 
behaviour. In this case, in spherical coordinates {G, and using the 'barred banal conformal representation' T = 4T, 
E + U = {E + U}/4, the time-independent Schrodinger equation is 

[Now E does not need redefining since the conformal term is zero because the configuration space dimension is 2.] While 
in plane polar coordinates {TZ, obtained by passing to stereographic coordinates on the sphere and then passing to 
the 'tilded banal conformal representation' T = T{1 + TZ^}^ and E + U = {E + U}/{1 + TZ^}^, the time-independent 
Schrodinger equation is 

' 'KS + ^SI = m) + U(R))». (13) 



TZdTZ { dTZj TZ^ 9$ , 

2.3 Inner products for scalefree triangleland in various coordinate systems 

The Jacobians in question are to be read off the metrics in Sec 1.3.2. In the barred banal conformal representation, 
this gives in spherical coordinates the well-known sin0 factor, in stereographic coordinates it gives TZ/{1 + TZ^}^ (which 
combines the usual TZ of plane polar coordinates with the requisite conformal factors), and, if one uses Ii (or I2 or ii or 
L2) instead of TZ, it gives but a constant factor. As there is no nontrivial weight in this representation, the above are the 
entirety of the extra factors inside the inner product (up to proportion, which one would then fix by normalization). In 
the tilded banal conformal representation, the Jacobians are the usual TZ in plane polar coordinates(7?., $), sec^^sinO in 
spherical coordinates (which combines the usual sinO with the requisite conformal factors), if one uses Ii, it gives 1/{I— Ii}'^ 
up to proportionality, and, if one uses I2, it gives l/l2- However this representation also has a nontrivial weight [c.f. (O], 
so that one ends up with precisely the same inner product as for the barred banal conformal representation (as should be 
the case, for in configuration space dimension 2, the wavefunctions themselves do not scale). 

2.4 Separability for ^-independent potentials 

Then the above time-independent Schrodinger equations are separable under the separation ansatze 

*(7^,$) = C(7^)?7(^') or *(e,$) = C(e)r/($) , (14) 

and in each case one obtains simple harmonic motion solved by 

77 = exp(±zj$) (15) 

for j an integer. This is a relative angular momentum quantum number corresponding to being classically conserved, 
in analogy with how there is an angular momentum quantum number m in the central-potential case of ordinary QM 
corresponding to the angular momentum being classically conserved. The accompanying separated-out equation is, in 
the tilded banal conformal representation in {TZ, <&) coordinates, the radial equation 

TZ^Cnn + TZCn - {2TZ^{\/{TZ) - E{TZ)}/h^ + j'}C = , (16) 
or, in the barred conformal representation in (O, <&) coordinates, the azimuthal equation 

{sine}-Hsiii0ee}e - {2V(e) - E}/h^ + f{sme}-^}^ = . (17) 

2.5 The special harmonic oscillator quantum problem 

The problem is, in the barred banal conformal representation in spherical coordinates, p2p with an harmonic oscillator 
type potential (1.59) inserted in, which I rearrange to the dimensionless form 

sine— [ + -^— + {£-A- Bcose}^ = (18) 



sine del de sin^e 
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for 

A = 2A/h^ , B = 2B/h^ , I = 2£/h^ . (19) 

This then separates into (fT5|) and 

{sine}-HsineCe}e - {A-£ + Scos(e) + j^{sine}-^}e = . (20) 

Alternatively, in the tilded banal conformal representation in [TZ, $) coordinates, it is (fT3|) with (1.58) inserted in, 
which 1 rearrange to the dimcnsionlcss form 

ndn\ dn] n-^d<^^- Xii + w^Y {i + n^y^] ' ^ ' 



for 

This then separates into p5|) and 



JC, = K,/h^ and £ = 2E/r . (22) 



«'^«»+«««-{Trw'(TW+-T = °- '''' 

Note that these equations are self-dual in the sense of so that, again, direct study of only one of the two asymptotic 
regimes is necessary and then the other can be read off by mere substitution. 

2.6 Useful mathematical analogue 

It is next convenient to recollect my observation [8] that scalefree triangleland's special triple harmonic oscillator like 
potential has the same mathematical form as the fairly well-known problem of the linear rigid rotor in a background 
homogeneous electric field in the symmetry-adapted basis for that problem. This then at the quantum level amounts to 
the present mathematics being analogous to that of the Stark effect for the linear rigid rotor, which is well-documented 
(see e.g. and Sees 2.5, 2.7, 2.8, 3). 

2.7 Exact solution for the 'very special multiple harmonic oscillator' 

In this case the analogy is with the linear rigid rotor itself [HI [52] , and the spherical representation's mathematics takes 
a familiar form, (1^0]) now being the associated Legendre equation (|75)) . 

One then has a quantum number J g No analogous to the total angular momentum quantum number 1 of the linear 
rigid rotor (the analogy is that J is a total angular momentum in a 3-d space, but that 3-d space is not ordinary 3-d space 
but rather the scaled triangleland configuration space, c.f App. l.C). This obeys 



and hence 

in the spherical representation, or 



J{J + 1} =£'-y^ = 2{E-A}/fi2 (24) 
E = ft2j{J + l}/2 + A (25) 



E = 2ft2j{j + 1} + {K^ + X2}/4 (26) 

in the planar representation. One also has a quantum number j G Z such that |j| < J; this is analogous to the magnetic 
quantum number of the linear rigid rotor, and here has the interpretation of being the 2-d relative angular momentum 
(which is the Z-component of the above 'total angular momentum'). 

Moreover, if E, i^i, K2 are to take the interpretation of being fixed, then there will either be 1 or such J - a closed 
universe type truncation of the number of allowed states. From (|26|) then, for there to be any chance of solutions one 
needs E > {Ki + K2}/'i, so for harmonic oscillator type models, E > is indispensable. If there is a J and it is not 
zero, there are various degenerate solutions corresponding to different values of j. These correspond to states of different 
relative angular momentum between the particle 2, 3 subsystem and the particle 1 subsystem. 

The wavefunctions are of the form 

*(e, $) cx Pj(cose)exp(ij$) . (27) 
Thus, in the planar representation, the solution is 



l-7^■ 



M'Jj(7^, $) « Pj ) exp(*j$) , (28) 

and, in terms of the physically and visually useful partial moments of inertia-relative angle variables, it is 
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t2, $) « Pj (^^^ exp{ii<f) = Pj f exp(zj$) = Pj f ^^V^ ) exp(jj$) . (29) 

While, in terms of the original variables of the problem, 

^liiL^,L2) « P\ f|4fe-K ) exp ( ^j arecos ( ) ) . (30) 



VM.llP + IU2lP 

[Figure 1. For the first three values of J (0, 1 and 2), I provide 

a) the azimuthal probabihty density function on the sphere, which is standard. 

b) The radial probability density function in the {Ti, <E>) plane Unlike the usual situation with the radial probability density function in the 
atom, the inner peaks are the taller ones (in the atom it's more probable e.g. that the 2s electron is 'outside' the Is one). This is due to the 
unusual inner product of our planar problem. 

c) These plots can then be interpreted in terms of straightforward relational variables such as (.1,12, II or I2. I plot in terms of Ii. Note the 
reflection symmetry about 1/2 and also that plots in terms of I2 coincide with those in terms of Ii and so I do not provide them. Prom the 
various above plots, one can infer on which of the sorts of triangles defined in Fig I.5d) the various wavefunctions peak. The cases depicted 
exhibit one of the following patterns of likelihood for configurations (triangles): a Jacobi-regular peak (if |j| = J), a Jacobi-tall and a Jacobi-flat 
peak with an interposed Jacobi-regular node (if |j| = J - 1), a Jacobi-tall peak, a Jacobi-regular peak and a Jacobi-flat peak with two interposed 
nodes of more moderate tallness and flatness (|j| = J — 2). Note that the higher J is, the more pronounced the tallness and flatness involved 
is. Next, note that for j = all the wavefunctions are surfaces of revolution. Higher values of j have sinusoidal dependence on <&. Moreover 
the directions picked out by this have physical meaning. |j| = 1 has what would usually be p^' and py (or dxz and dyz) directionality in space, 
which in our problem signifles near-coUinear peaking and near-isosceles peaking in the conflguration space of triangles. While, |j| = 2 picks out 
both of the above equally (dxy) or avoids both equally {d^2_y2).] 

Note 1) paralleling the classical working, the very special solution is unconditionally self-dual under the duality map. 
2) while this case is simple and not general, its exact solution serves as something about which one can conduct a more 
general perturbative treatment (SSec 2.10). 

2.8 Small asymptotics solutions for the 'special multiple harmonic oscillator' 

I work in the tilded Q-representation. In the first small approximation, 2{E+ U} = Qq (a constant evaluated in Sec 1.4.9). 
Thus, in terms of Qq = Qq/H^, the time-independent Schrodinger equation is 

which is a familiar problem, separating into simple harmonic motion p5p and a radial equation (jl6p . which is now the 
Bessel equation of order j ([55]) in the rescaled variable S = ^/QqTZ. Thus there is just the one quantum number j G Z 
(relative angular momentum). The wavefunctions are then 

^'j(7e,$) cx Jj(yQ;;7^)exp(^j$) . (32) 

So, in spherical variables, 

and in terms of the partial moments of inertia-relative angle variables, it is 



*j(e, $) cc Jj(v^tan|)exp(ij$) , (33) 



*j(Ii,l2,$) cx Jj(VQoIi/l2)exp(zj$) = Jj(VQoIi/{I-Ii})exp(zj$) = MV Qo{l - h} /h)exp{ti<^) . (34) 
While, in terms of the original variables of the problem, 

*j(tl,i2) « -^j ( 



exp(^zjarccos(^-^jj . (35) 

I sketch the probability density function in various of these variables in Fig 2. Finally note that replacing the specific 
constant Qq by some other constant Q, this also provides the exact solution to the constant potential problem. 

In the second approximation, 2{E -I- U} = Qo ^ Q2T^, where Q2 is another constant evaluated in Sec 1.4.9, and which 
I consider here to be strictly positive. In this case the time-independent Schrodinger equation is 

(where I use Q2 = Q2/fi^), which is also a familiar equation, straightforwardly mapping to the 2-d isotropic harmonic 
oscillator [5011541155] by the correspondence in App B. It separates into simple harmonic motion (|15p and a radial equation 
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that this time can be mapped to the associated Laguerre equation Thus this problem's solution is as follows. 

In addition to the relative angular momentum quantum number j G Z, there is a 'radial', 'principal' or 'node-counting' 
quantum number R g No such that 

Qo/2^V^= l.il + 2R+l (37) 
holds. In terms of the original quantities of the problem, this corresponds to 



E = K2/2 + hy/4E + Ki 
or, in the language of the spherical problem 



3X2 



2R+ 1} 



A + B + hVE - A-2B{ |j| + 2R + 1} 



Then 



E > K2/2 + Hy/AE + Ki - 3K2 ( >0) 



(38) 



(39) 



(40) 



is necessary for there to be any solutions at all. While, the non-standard interpretation that all of E,Ki,K2 are fixed 
will restrict the number of simultaneously-relevant solutions in my closed-universe context; e.g. rational-irrational incom- 
patibility can be used to construct cases with no solutions, but there are also cases for which there remain interesting 
degenerate possibilities such as (R, j) = (0, ± 2) and (1, 0). 
The wavefunctions are 

*Rj(7^,$) cx 7^ljlcxp(/&7^V2)41(v^7^2)exp(y$) 



or, in spherical variables, 

^'Rj(e,$) cx tanljl|exp(Q2tan2|/2)Ll:['(v/Qltan2|)exp(ij$) , 
while in terms of the partial moments of inertia-relative angle variables, it is 



(41) 
(42) 



«'rj(Ii,l2,$) oc exp 



2 I, 



Ii 



I-Ii 



|j|/2 



exp 



Ii 

2 I-Ii 



r I- 



■I-Il 



cxp(ij$) = 



I-I2 



|j|/2 



S2^^ exp(ij$) = 

\/&I-l2 



exp 



-I-I2 



While, in terms of the original coordinates for the problem. 



|j| 



exp 



7-lj 



exp I] arccos 



Ll • k2 



exp(ij$) 
(43) 

(44) 



That a 2-d isotropic harmonic oscillator with eigenvalues (P71) resides within the scalefree triangleland multiple har- 
monic oscillator like problem as a limit problem has counterpart in the literature on the linear rotor (|56j. see also the 
figure in [57j). 

[Figure 2 a) Sketches of probability density functions for the first approximate solutions in terms of TZ for j = and j = 1 (which is qualitatively 
similar to j = 2). For Qo = 1, what the relevant region contains is but the first bit of ascension to the first of an infinity of peaks (first two 
subfigures). For Qo large various whole peaks can lie in the relevant regime; the two plots I provide are for Qo = 10* with j = and j = 1. 

b) Sketch of the (unnormalized) probability density functions for the first approximate solutions in terms of Ii for j = and j = 1, for Qo = 1 
and 10'*. Note that for the first approximate solution in terms of I2, all one need do is reflect in 1/2 (and the relevant region maps to the new 
relevant region and so is the same portion as before). Thus I do not provide separate figures for these but rather provide the whole range of 
values of Ii. This refiection trick works throughout this Section. 

c) The second approximate solutions in terms of TZ for (R, j) = (0, 0), (0, ±1) or (0, ±2) and (1, 0) for Q2 = 1 and then for Q2 = 100 in h = 1 
units. 

d) The second approximate solution in terms of Ii for these same six cases.] 



Note 1) As regards which inner product to use, in the first approximation, one can use the naive TZ in place of TZ/{1 + TZ^}'^ 
at the cost of an error of up to 2% over TZ € [0,0.1]. This however swamps the improvement in passing from first to 
second approximation; one should use at least TZ{1 — 2TZ^} in the latter case so as to not compromise accuracy. 

2) The unusual inner product in use here does cause these probability density function differ somewhat from the usual 
ones for a 2-d constant potential or 2-d isotropic harmonic oscillator: the peaks are somewhat shifted inwards, and the 
inner peaks gain in relative importance; these are signs of greater confinement. 

3) How can few-peak functions be second approximations when infinite-peak functions are first approximations? The 
trick is that the few-peak functions have an extra node-counting quantum number R, so that the first approximation can 
be seen as a superposition of many different values of R each with its own Q2 as dictated by ([57)l . and thus build up a 
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multiplicity of peaks. For the values considered, there is better than 1% accuracy between the 2 approximations up to 
some value of TZ of order of magnitude 0.01 to 0.1. 

2.9 Large asymptotics solutions for special harmonic oscillator and other problems 

I now apply duality to write down first and second large approximation solutions. The first approximation gives 
^'j(7^,$) cx Jj(v/QI/7^)exp(y$) or «'j(e, $) cx Jj(v/Qlcotf )exp(iij$) . 



' 2 . 

Thus, in terms of the physically and visually useful partial moments of inertia-relative angle variables, it is 



*j(Ii,l2,$) = Jj(VS4l2/Ii)exp(zj$) = Jj(v/Q4{I-Ii}/Ii)exp(zj$) = Jj (v/Q4l2/{I - l2})exp(zj$) . 
While, in terms of the original coordinates of the problem, 

ti • t2 



*jUl,A2) OC Jj 





L2 






^ exp ^ 









|'-l||||'-2| 



(45) 



[Figure 3 Using the same 'radial' inner product as in figure 1, I obtain tlie following probability density functions in terms of TZ. [In terms 
of Ii , these have the same form as the corresponding small solution's probability density function in terms of I2 and vice versa, so I do not 
provide any new probability density functions as functions of Ii, I2.] 

a) For the first approximation, I provide an archetypal probability density function for Q4 = 1 and 10*. 

b) For the second approximation, I provide the probability density function's for (R, j) = (0, 0), (0, 1) and (1, 0) [{0, 2) is qualitatively similar 
to (0, 1)] for Qe = l and then 10"'.] 

As happens classically, the above analysis also holds for Vj^, q) with constant of proportionality A^, one gets exactly 
the same large-asymptotics analysis as here, with qq = 2E — Aq, (72 = 4E — {4 + a}Aa. Thus the conditional duality map 
has put the curious large asymptotic region in firmly understood terms. 

2.10 Special triple harmonic oscillator problem treated perturbatively for small B 

Let us treat the special triple harmonic oscillator problem as a perturbation about the simple 'very special triple harmonic 
oscillator'. For a perturbation Ti' to one's rescaled Hamiltonian [rescaled as the other calligraphic quantities in (|19p are] 

the first few objects of perturbation theory as used in this paper are as follows [52]. For £j ~ nondegenerate, — 

{J\n'\3), = - Ek#,j ^^tS!^^ at first order and = - ^^^^j ^^^^Iff ^^'^^^^ ^J°^ degenerate, 

one needs to solve at first order X]s(J' Jl^'l J' s)ai^^ = £^^''a\^\ and, at second order, J2k^] J2s ^^'^I^ |J,s)(j,s|'H |K,k) ^(2) _ 

For us, the = |J,j) are ^7}i and H' = B'X for X the 'Legendre variable' related to 6 by X = cosO. Then the 

key integral underlying time-independent perturbation theory is (^fj/j/ |V'|5'jj) (here there's no subtlety in the style of Sec 
2.3 with the inner product as the banal conformal transformation here has a merely constant conformal factor and so 
cancels out upon performing normalization). The above has as its nontrivial factor J -^ Py{X)XP'^{X)dX , but there is 

a recurrence relation ([78]) enabling XP\{X) to be turned into a linear combination of Pj„(Ar), whereupon orthonormality 
of the associated Legendre functions ([77]) can be applied to evaluate it. This calculation parallels that in the derivation 
of selection rules for electric dipole transitions [5T1[5H]. Then (J, j'|KAr|J, j) = since the recurrence relation sends XPj^ 

to a sum of Pj] ' for j' ^ j, so each contribution to the integral vanishes by orthogonality. While, for (J', j'jPArjJ, j), one 
similarly needs J' = J ± 1 and j' = j to avoid it vanishing by orthogonality. So two cases survive this 'selection rule'. The 
first is, by direct computation, 

while the second is 



(J - 1, jl^XIJ, j) = (J, jl^XIJ - 1, j) = bJ ^.jJ^I^^j^,^ (47) 



by using (|46)) with J - 1 in place of J (paralleUing e.g. [50]'). 
The eigenspectrum is thus 
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so 

ft2j{J + l}{2J-l}{2J + 3} 
For J = 0, one needs a separate calculation, which gives 



Eo,o=4A+ — +0(^4) . (50) 

This calculation can then be checked against its rotor counterpart (originally in [59 and which can also found in 
e.g. |49[ 150] ]. The corresponding eigenfunctions can be looked up (e.g. [60]) and reinterpreted in terms of the original 
problem's mechanical variables. 

2.11 Placing a closed-universe interpretation on the perturbed problem 



N.B. the above are universes not modes within a particular universe - one has just the one value of E, Euniv. Sometimes 
|19| . this corresponds to no allowed J, j, sometimes to one and at least sometimes to more than one. For example, 
£^0,0 = £'1,0 for B — y^l5/2h^, which is perturbatively acceptable provided that A >> B. From this, one can extract 
partition functions using [61j to take into account that the energies are only known perturbatively. Hence one can extract 
a classical notion of entropy and hence of information. Also, from knowing the wavefunctions, one can construct mixed 
states and then the quantum-mechanical von Neumann information corresponding to these. 

But for records-theoretic purposes, it is information content of subsystems and relative information between subsystems 
that look to be more significant quantities. These are obtained from solving the subsystem quantum problems (subject to 
global restrictions such that the subsystem energies add up to a fixed energy of the universe [5]) and then proceeding to 
compute notions of information similar to the above (but now including relative information between pairs of subsystems 
that is a quantity of higher relevance as regards records-theoretic approaches). This would require considerable further 
study. 

2.12 Exploiting further correspondences with the rotor problems 

The following are available in the rotor literature (in the spherical presentation) and could thus be straightforwardly 
transcribed to the various presentations for this paper's scalefree triangleland problem. Higher order corrections in B 
are in the literature for the rotor [62l [63] and so can be transcribed into the relational context [e.g. I was able to write 
0{B'^) and not 0{B^) in l|49p due to this]. Alternative variational methods (using the Hellmann-Feynman and hypervirial 
theorems) appear in |63j : these cover higher order terms too, and can be used to show generally that only even powers of 
B occur. The calculation for the large B regime has both been done |SD] and matched to small B regime calculations. 

Numerical evaluation of eigenvalues was been done by Lamb's [M] continued fraction method [511[nSl[5B] (and otherwise, 
e.g. |57l [63]). Myself, I just simply used Maple's [53] rkf45 solver alongside the iteration of the method in e.g. [55] to 
locate the eigenvalues, which gives reasonable agreement with the formula at the end of SSec 2.10. 

The rotor literature also indicates how the C perturbation can be transformed away with a new rotated choice of 
coordinates in which the maths is again that of a B type perturbation (Sec 3 - the quantum application of the trick in Sec 
1.5). While in the laboratory with a rotor one could choose one's axial 'z' direction to be in the most convenient direction, 
there are various Problem of Time strategy modelling reasons not just to stay in these coordinates in our triangleland 
problem (see the next Sec). 



3 General triple harmonic oscillator problem for scalefree triangleland 

In the spherical coordinate representation, this takes the form 
where A, B, C are given by and 

C ^ 2C/h^ . (52) 

This is harder because the angle-dependence of the potential results in nonseparability in these natural coordinates, which 
makes for a useful model of the semiclassical approach to the Problem of Time [47 . 



3.1 Solution via use rotation/normal modes/adapted basis 

As in [8], start again using rotated/normal coordinates at the classical level, or switch to such coordinates at the differential 
equation solving stage, amounting to a choice of basis in which the perturbation is in the (new) axial 'zn' direction. While 
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it can be viewed as before in the new rotated/normal coordinates, nevertheless studying the original coordinates' $- 
dependent V term remains of interest as it may well be appropriate for the original coordinates to have mechanical 
attributes or the heavy-light subsystem distinction underlying the semiclassical approach!! Moreover in that case, being 
able to proceed further in the rotated/normal coordinates can serve as a check on "standard" procedures in the original 
coordinates (along the lines suggested in |37j). e.g. as a test of the validity and accuracy of assumptions and approximations 
made in the semiclassical approach. This section serves to begin to set up such a check-point. Another reason why one 
might not adopt the normal coordinates as the physically significant ones is in order to use J as a trackable non-conserved 
quantity as regards investigating the semblance of dynamics in a fundamentally timeless records theory approach. 

There is now an issue in the projection that there is an additional factor due to the change in area in moving between 
each patch of sphere and each corresponding patch of plane. Namely, the probability density function on the sphere is 
sin8|5'(0, <I>)p of which the sin0 pertains to the sphere itself, while the probability density function on the stereographic 
plane is 7?./{l + TZ^}\'i'{TZ,^)\^ of which the Tl/{1 + TZ^} pertains to the stereographic plane itself. The very special 
solution now sufhces as an illustration. The analytic form of its probability density function is 



PDF(eN, *n) oc sineN{Pi(coseN)}^ 



(53) 



and so 



and so 



-nr^T^rr^ ^\ ■ ^ I f ScosO -f CsinGcos* \ 
PDF(9, $) cx sm6 Fj ( — j 



PDF(7^, $) cx 



n 



s{l-7^2} + 2c•7^cos$^ 



1 + 7^2 

Then, in terms of the partial barycentric moments of inertia, 



SN{l + 7^2} J 



(54) 
(55) 



probability density function(Ii, $) cx < Pj 



B{1 - 2Ii} + 2CVli{I-Ii}cos$\ 



(56) 



and 



(57) 



For (J, j) = (0, 0), these are the same as for the C — case, so there is no need to provide a new plot. However for higher 
values of (J, j), they are distinct (Fig 4). 



[Figure 4 For B = 1 with C = 0.1, then 1 and then 10, I give the following triples of plots. PDF{7^, *) for (J, j) = (1, 0). PDF{n, *) for 
(J, j) = (2, 0). PDF{Ii, <!>) for (J, j) = (1, 0). PDF(Ii, <!>) for (J, j) = (2, 0), the first of which has an shallower ring hidden inside the visible 
outer ring.] 



4 Conclusion 

4.1 Results Summary 

Study of relational particle mechanics is motivated by the absolute versus relative motion debate and the analogy between 
relational particle mechanics and the canonical formulation of General Relativity. Specific quantum similarity relational 
particle mechanics of N = n -I- 1 particles in 2-d can be constructed due to knowledge in this case that the classical 
configuration space is CP"^^ [7l[TTl[66]. Likewise, specific similarity relational particle mechanics of N particles in 1-d 
can be constructed due to the knowledge that in this case the classical configuration space is S"^^. This paper provides a 
quantum study of similarity relational particle mechanics with harmonic oscillator like potentials, for which I use exact, 
asymptotic, perturbative and numerical methods. Throughout, a mathematical analogy with the linear rigid rotor in 
a background electric field is useful; I then transcribe this from spherical/stereographic plane terms to be in terms of 
(partial barycentric moments of inertia, relative angle) variables and the original similarity relational particle mechanics 
problem's mass-weighted relative Jacobi variables. In particular, 

1) in the spherical representation there is a constant-potential subcase soluble thereupon in terms of spherical harmonics. 

2) I then consider relative angle $ independent potentials, corresponding to the existence of a relative angular momentum 
type conserved quantity whereby the classical and quantum theory is simplified. This permits solution of the large 
and small stereographic coordinate regimes in terms of Bessel functions, and, more accurately, in terms of Laguerre 

^In the laboratory, one might likewise not pick the normal coordinates of the rotor— electric field if there is e.g. also a magnetic field that 
picks out a different direction. 
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polynomials; the large regime is, moreover, universal rather than dependent on the choice of harmonic oscillator like 
potentials. 

3) There are various theoretical reasons to wish to remove the $-indepenent restriction - dynamical and quantum me- 
chanical nontriviality and genericity, as well as toy-modelling semiclassical and records theory approaches to the Problem 
of Time in Quantum Gravity. Thus I also treat ^-dependent potentials by a coordinate rotation/normal modes/adapted 
basis construction. 

4) In each case, I then interpret these solutions in terms of the underlying mechanical variables, including investigation 
of which triangles formed by the particles are more and less quantum-mechanically probable in a given state. 

4.2 Further Extensions 

The present paper is furthermore useful in that many methods used in it can furthermore be used toward solving other 
concrete relational particle mechanics examples. Following the success of identifying the scalefree triangleland multiple 
harmonic oscillator like potential problem with the linear rigid rotor Stark effect, surveying the Quantum Chemistry 
literature for analogues of some of the below may be useful. 

Scalefree 4-stop metroland also has the configuration space S^. 

Then the conformal-ordered time-independent Schrodinger equation for multiple harmonic oscillator potentials is 

1 d ( dW 1 1 (9^vl> 

iiine—} + —^--r=A + Bcos{2Q)+Csm^ecos{2<^>) , (58) 



sine de\ dej sin^e 
or, in Legendre variables, 

dxV^- ^'^dx\ + T^^d^ = ^ + ^'^^^ ~ ^^="^(2$)} Ccos(2$) (59) 

where A = 2{A- E}, B = 2B/fi\ C = 2C, andV = A- B. 

Then for S = C = 0, one obtains the Legendre equation. Then for the B perturbation, double use of the recurrence 
relation ((77)FI gives rise to a Aj = 0, AJ = 0,±2 'selection rule'. While, the C perturbation has 

j\ Pj=^^(X){l - X'^}Pl{X)dX contributions, for which double use of the recurrence relation ([78|) gives in this case the 
selection rule Aj = ±2, AJ = 0, ±2. Note that normal mode trick doesn't work for this problem, so here one has to study 
the C perturbation as well as the B perturbation. 

Higher N-stop metrolands 

Here the conformal-ordered time-independent Schrodinger equation is, in terms of ultraspherical angles, 

where £ — 2E/?i^, JCp — Kp/h?, rip is the unit vector of the embedding Euclidean configuration space 7?.(N, 1) = R" and 
E has been redefined to incorporate the conformal term since (hypcr)spheres are of constant Ricci scalar. 

There is also a highly special constant potential case within the multiple harmonic oscillator like potentials. This is 
now a more complicated sequence of associated Gegenbauer equations as explained in Appendices A and B; these are 
nevertheless also fairly standard and well-documented [671 168] V One can then study perturbations about this. Then, if 
the associated quantum number is not zero, one does not get the Gegenbauer pairings or the right weights straight away 
(see Appendix B). Computation of first order perturbations in this case requires the Gegenbauer parameter converting 
recurrence relation (jSip as well as the polynomial order reducing recurrence relation (|80p. making the calculation somewhat 
more complicated in this case. 

Yet further extensions 

As regards (N > 3)-a-gonland's 'genuine' CP*^ mathematics (rather than CP""^ mathematics that is re-expressible as 
mathematics) is required, placing this beyond the scope of the present paper. 

An alternative type of extension is to keep a given geometry but consider a variety of other potentials, e.g. Coulomb-like 
potentials or a mixture of Coulomb-like and harmonic oscillator like potentials. 



'An alternative computational scheme is to proceed by the formulae for integrals of products of three spherical harmonics in e.g. ([58]). 
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4.3 Applications to classical and quantum geometrodynamics 

The semiclassical approach to the Problem of Time in Quantum Gravity requires nonseparabihty so that the crucial cross- 
term in ([7]) exists. To set up such a model, one could have e.g. two heavy (H) particles and one light (L) one, leading to 
one H relative Jacobi separation and one L one. The nonseparabihty requirement would mean requiring the interpretation 
that the physics picks out 8 and $ variables unaligned with the simplifying normal modes ones (in terms of which there 
is separability), and that the potential be ^-dependent. Another way of setting up a semiclassical approach model using 
the material in this paper is if O and $ are considered to be H and L respectively (the opposite assignment is impossible 
as sin^O cannot be >> 1). The above two situations would ba a significant improvement on the example in [37] because 
that has no nontrivial linear constraints; it remains to be seen how far such a calculation could be taken analytically (and, 
if needs be, numerically and/or subject to further approximations). While, the [37] example's ability to be solved by 
usually-unavailable means external to the semiclassical approach is retained by the above examples - it is the work in the 
present paper. Thereby one can assess whether the semiclassical ansatze and associated (and any other non-associated 
but unavoidable) approximations are sensible for these new semiclassical approach models. To have a shape-scale aligned 
H-L split, which in some ways more closely parallels GR Cosmology (H scalefactor versus L inhomogeneities and/or 
anisotropics), one needs to study the Euclidean relational particle mechanics counterpart of this paper |17| . This is also 
the case if one wishes to obtain a dilational (York- like) internal time model [HI • 

As regards building concrete relational particle mechanics examples of the records theory approach to Quantum 
Gravity, a notion of distance on configuration space readily follows [40] from the metrics in As regards computing 
a notion of information/negentropy, given an explicitly solved QM, one can (c.f. Sec 2.11) build a statistical mechanics 
from that and extract the negentropy/information. Moreover, this continues to be the case when the QM is only known 
perturbatively [61] (with correction terms up to the corresponding perturbative order [611), which is one reason for the 
relevance of the perturbative calculations in this paper. Further questions then are: does this paper's model have a notion 
of information storage? (In parallel with [39], might a heavy particle passing by ~ the signal - imprint separation/motion 
on the other two particles - the record? [39] proceeds to study this via path integrals and associated objects such as the 
influence functional; the extent to which these more specialized computations have been carried out for the analogous 
rotor is an interesting question, whose investigation might substantially cut down on how much work is needed to complete 
the parallel calculation to [39] for the present paper's models.) 

As regards semblance of dynamics emerging from timeless records theory, the relative angular momentum that 
was a conserved quantity for B = becomes a changing quantity for C 7^ 0, so tracking and explaining that may be of 
significance. I can track this classically by e.g. computing it at each stage in the rkf45 routine. An interesting question then 
is whether and how this could be tracked quantum-mechanically? As regards whether evidence can be found for/against 
Barbour's conjecture of time capsules [24', '3], the probability density functions plotted in this paper are the right output 
for addressing that. The fairly standard maths I obtain (at least in my simple specific example of harmonic oscillator like 
potentials) suggests that, if time capsules do occur, they ought to be findable also within standard QM (for perturbed 
linear rigid rotors). However, further scaled triangleland or trihaedronland (— 3-haedronland) work would be necessary 
as regards more specific conjectures Barbour makes about time capsules that were specifically made about 3-particle 
Euclidean relational particle mechanics and whether its triple collision or uniform (i.e. equilateral) configurations play 
highly dominant roles (where the probability density function 'mist' might be highly concentrated). For the moment, my 
simple similarity relational particle mechanics model would not seem to exhibit very heavy peaking around its equilateral 
configuration. 

This paper's model is also conceivably an interesting one from the perspective of histories theory and as regards the 
problem of finding (partial) observables for Quantum Gravity. 

On the whole, this paper downplays the suggestion [TH (3^ that (Barbour's) relationalism requires radically different 
QM theory, though one would need to check further examples (including more complicated ones) to be more sure of 
this conclusion (some of the further examples above can be motivated on such grounds). Some closed- universe and finite- 
universe effects are, however, manifest. As regards the issue '19Jf5| of whether stationary quantum universes have single or 
multiple states, I comment that this paper's models do exhibit some degenerate states (both among simple exact solutions 
and perturbatively to second order); however one needs a much more extensive study of relational particle mechanics 
model universes before one can begin to say whether these are, however, non-generic. There is also a certain amount 
of tension between results obtained by reduced quantization as in this paper and by Dirac-like quantization as in [5], 
though I postpone discussion of this to [17 (one of the problems being that my method of Dirac-like quantization in [S] 
for Euclidean relational particle mechanics does not directly extend to similarity relational particle mechanics, so that 
checks between the two methods and lessons drawn from them are best left to the Euclidean relational particle mechanics 
arena) . 

Treating cases with more than three degrees of freedom will be necessary to get a grip of some aspects. This is the 
case firstly for investigation whether there are kinetic effects of the kind that Barbour conjectures [24. ^3] as regards a 
semblance of time arising from timeless configurations. All of the §"^^ being conformally flat, one would need to extend 
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to such as the non-conformally flat CP^ so that eS'ects that are irreducibly kinetic occur (rather than cases for which 
a kinetic and potential re-definition leaves one with a flat kinetic term in the end by passing the conformal factor into 
a re-defined potential). Secondly, for (N > 4)-stop metroland, the relative angle dependence comes from writing the 
potential in coordinates dictated by the kinetic term, and no longer in such a way that the standard rotation to normal 
coordinates removes this complication. A third such point, which constitutes an interesting further investigation in its 
own right, is that N-particle models also have a robustness application: is the QM of N - 1 particles stable to the 
inclusion of a further particle? This parallels the situation of whether Taub space is stable within the mixmaster solution 
in minisuperspace Quantum Cosmology (found to be unstable in |69|V Fourthly, upgrading to models with > 3 particles 
is likely to improve one's capacity to use particle clumps to model inhomogeneities. A final question that such models 
can be used to investigate is whether the conformal ordering succeeds in avoiding confiict with other necessary technical 
conditions. 
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Appendix A 

Following from App I. A. 2, the energy constraint gives, in conformal ordering, the time-independent Schrodinger equation 



A.l Scalefree N-stop metroland 



Here, by p 1269-70 of JX)j, an appropriate finite subalgebra acting on the corresponding cotangent space is then SO(nd)(S)R"'^, 
where (§) stands for semidirect product. This can be considered to be generated by angular momenta Jj and coordinates 
Ui such that = 1- 

Now, the Laplacian corresponding to line element (1.20) is 

(62) 

This situation applies to preshape space and for shape space [66l HI] ■ Hence QM on preshape space for N particles in 
dimension d has the time-independent Schrodinger equation 

(where the energy has been displaced by a constant term from the constant-curvature contribution to the conformal order- 
ing). The d = 1 case of this is, additionally, the time-independent Schrodinger equation for scalefree N-stop metroland, i.e. 
eq. pO| . These Hamiltonians involve suitable quantum operators (see e.g. p 160 of [71] for an account of the properties 
of the constituent Laplacian operator on S"~^, see also [72j[73]). 

For constant potential, this time-independent Schrodinger equation is an equation of form D^'i> = A*!*. Then the 
separation ansatz ^' = 11^=1 ■0p(0p) yields the simple harmonic motion equation for Qn-i and n - 2 equations of form 

n v'2 id^V^n-p ; ..^ dljjn-p , . r. , ol / jp-2{jp-2 +P- 3} 

{1 - Xn-p}-r^ {P - ^}^'^-p-n? + Jp-l{jp-l + P - 2}V'n-p z -2 i^n-p = (64) 

under the transformations — cosGp, p = 1 to n - 2. These are associated Gegenbauer equations (|82[) with parameter 
Xp = {p— 2}/2, where the integers jp_i, jp_2 G No are picked out as eigenvalues, so that the {n - p}th equation is solved 
by C-''"^ (cos0„_p; {p — 2}/2). Then one gets a sequence of integer quantum numbers beginning with the familiar |ji I < jo- 
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A. 2 Scalefree N-a-gonland 

In this case, as regards kinematic quantization for scalefree N-a-gonland, use that S(N,2) = CP""^ ==SU(n)/U(n- 1), 
which is of the general form Q ^ G/H considered in [70] . Thus a suitable finite algebra acting on the corresponding 
cotangent space is SU(n)@R^". 

The Fubini-Study Laplacian for S(N, 2) is then 



= ' " " -' <^ — — <^ '■'■P=' - i^P? -I- VPTin—- V + — — <^ ""P^" - <^ — + 1 

^;=^7^p \97^p \{l + ||7^||2}2n-3l'' ^ f ^ dQp\{i + \\n'V''-' Wp ' i ae^ 

■ (65) 

so that the conformal-ordered time-independent Schrodinger equation for the 2-d shape space of N particles is pT|) . 
A. 3 Scalefree triangleland case 

In the tilded banal conformal representation, the (3, 2) case's Laplacian in the flat coordinates {TZ, <&) takes the familiar 
form 



While, in the barred banal conformal representation in spherical coordinates (0, <&), it takes the also-familiar form 

^ =^7i7^isine— ^ + -^— 5 . 67 
sm9 dQ [ 96 J sm^^S 9$^ 

While I provide a parallel scheme for general CP*" above, it is the CP^ = version actually used in this paper's calculations 
for which I currently have guarantees of good behaviour as an operator. 

Appendix B: Special functions and mappings of ordinary differential equations 

The Bessel equation of order p, 

v'^Wvv + vwy + {t;2 — P^}w — , (68) 

is solved by the Bessel functions. I denote Bessel functions of the first kind by Jp{v). 
The associated Laguerre equation, 

xVxx + {a + l~ xjy^ + ny = , (69) 

is solved by the associated Laguerre polynomials ^"(x) (and unbounded second solutions). The 2-d quantum isotropic 
harmonic oscillator's radial equation for a particle of mass fi and classical oscillator frequency w, 

maps to the associated Laguerre equation under the asymptotically-motivated transformations 

R = {hx/niu}^ e'^'/^yix) , x = fnur^/h (71) 

and so is solved by 

corresponding to the discrete energies E = {|m| + 2r + l}hLO for radial quantum number r G Nq |54l 155] . 
The associated Legendre equation 

{1 - X^jYxx - 2XYx + {J{J + 1} - j'{l - X^}-'}Y = (73) 

(which is equivalent to the equation 

sin-ie{sineye}e + { J{ J + 1} - j^sin-^ejy = Q (74) 

under the transformation X = cos8), is solved by the associated Legendre functions pj^'(X) (and unbounded second 
solutions), for J e No, j G |j| < J. We use the standard convention that 

P{(X) = {-l}'(l-.Y^)ig|T{5^;Ta|lf-'^'-l}'} ■ W 
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whereupon 



+ (76) 



\\ 2 {J + |j|)! , 

is a complete set of orthonormal functions for X E [-1, 1]. We also require the recurrence relations [68l \67\ 



XP, {X) ^^^^ , (77) 



Vl - X^P\-^ = . (78) 

The Gegenbauer equation 

{1 - X'^}Yxx - {2A + l}XYx + J{J + 2A}r = (79) 

is solved boundedly by the Gegenbauer Polynomials Cj{X\X). Normalization for these is provided in e.g. [BTlEH]; the 
weight function is {1 — 2 between equal-A Gegenbauer polynomials. These furthermore obey the recurrence relations 

[571 EH] 

^r(Y\\ {J + l}g.T+i(X; A) + {2A + J - l}g.,_i(X; A) 

XC,{X, A) = , (80) 

r ^^ A{Cj+i(X;A + l)-Cj_i(X;A + l)} 

cj+i(x,A) = jT}:ri ■ ^^^^ 

The associated Gegenbauer equation 

{1 - X^}Yxx - {2A + \}XYx + J{J + 2A}r - j{j + 2A - 1}{1 - X^y^Y = (82) 

is solved boundedly by the associated Gegenbauer functions Cj(X; A). These are re-expressible in terms of Gegenbauer 
polynomials via [72] 

C](X;A)«{l-X2}^Cj_j(X;A/2-l+j) . (83) 

With this conversion, the recurrence relations between Gegenbauer polynomials ([501 [ST]) turn out to suffice for this paper. 
For A = 1/2, ^ is the Legendre equation solved by Pj(X) = Cj(X;l/2), ^ becomes d??]), and ^ becomes the 
associated Legendre equation (|751) solved by Pp cx Cj(X; 1/2), by (|551) and [57] 

C,j(X; A) (X {X^ - \}-^--^p};l^{X) . (84) 
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